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A  Method  of  Solution  for  Painleve  Equations: 
Painleve  IV,  V 


by 

A.S.  Fokas,  U.  Mugan  and  M.J.  Ablowitz 

1.  INTRODUCTION, 

The  mathematical  and  physical  significance  of  the  S’x  Painleve 
transcendents,  PI-PVI,  has  been  well  established.  Their  mathematical 
importance  originates  from  the  following:  i)  P.  Painleve  [1]  and  B. 

Gambier  [2],  at  the  turn  of  the  century,  classified  all  equations  of  the 
form  qtt  *  F(qt,q,t)  where  F  is  rational  in  q t ,  algebraic  in  q  and 
locally  analytic  in  t,  which  have  the  Painleve  property,  i.e.  their 
solutions  are  free  from  movable  critical  points  [3].  Within  a  Mobius 
transformation,  they  found  fifty  such  equations;  these  equations  can 
either  be  integrated  in  terms  of  known  functions  or  can  be  reduced  to 
one  of  the  six  Painleve  ‘ranscendents .  ii)  R.  Fuchs  [4]  and  R.  Gamier 
[5]  considered  Painleve  equations  as  the  i somonodromic  conditions  for 
suitable  linear  systems  with  rational  coefficients  possessing  regular 
singular  points.  In  other  words,  Fuchs  and  Gamier  established  the  existence 
of  compatible  linear  systems  which  turn  out  to  be  the  analogue  for  Painleve 
equations  of  the  so  called  Lax  pairs  [6]  for  solvable  nonlinear  evolution  equa¬ 
tions.  The  condition  of  isospectral ity  is  now  replaced  by  isomonodromici ty .  How¬ 
ever,  apparently  the  above  authors  did  not  pose  the  question  of  using  these  Lax 
pairs  to  integrate  the  Painleve  transcendents .  iii)  Ablowitz,  Ramani  and  Segur 
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[7]  discovered  a  remarkable  connection  between  equations  with  the  Pain¬ 
leve  property  and  nonlinear  POE's  solvable  by  the  inverse  scattering 
transform  (1ST).  Invariant  solutions  of  such  PDE's  satisfy  equations 
with  the  Painleve  property.  (More  precisely,  associated  with  such  a  PDE 
consider  an  ODE  which  describes  those  solutions  of  the  PDE  which  remain 
invariant  under  the  action  of  some  Lie-point  group;  then  all  solutions 
of  this  ODE  which  can  be  obtained  via  the  1ST  of  the  corresponding  PDE 
have  the  Painleve  property,  see  also  [8]).  For  example,  proper  exact 
reductions  of  the  Korteweg-deVries  (KdV)  equation  lead  to  PI  and  PII 
[9];  PII  and  special  cases  of  Pill  and  PIV  can  be  obtained  from  the 
exact  similarity  reduction  of  the  modified  KdV,  sine-Gordon  and  the  non¬ 
linear  Schrodinger  equations,  respectively  [10];  special  cases  of  PVI 
can  be  obtained  from  exact  reductions  of  the  three-wave  resonant  inter¬ 
actions  [11]  and  from  the  Ernst  [12]  equation  [13].  iv)  H.  Flaschka 
and  A.  Newell  [14],  M.  Jimbo,  T.  Miwa  and  K.  Ueno  [15]  considered  Pain¬ 
leve  equations  as  i somonodromic  conditions  for  suitable  linear  systems 
possessing  both  regular  and  irregular  singular  points.  These  systems 
appear  more  suitable  than  the  linear  systems  introduced  in  [4]  and  [5] 
for  both  integrating  the  Painleve  transcendents,  as  well  as  for  studying 
their  asymptotic  behavior  [16].  Gamier  [17]  also  considered  linear  equa¬ 
tions  with  irregular  singular  points  in  connection  with  the  Painleve 
equations  but  apparently  did  not  pose  the  question  of  using  this  connec¬ 
tion  to  integrate  these  equations. 

The  physical  significance  of  the  Painleve  transcendents  follows  from 
their  applicability  to  a  wide  range  of  important  physical  problems. 
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Painleve  equations  appear  in  nonlinear  waves  (see  iii)  above)  in  quantum 
field  theory,  in  statistical  mechanics  [18]-[30],  etc. 

Associated  with  the  exact  integrabi 1 i ty  of  Painleve  equations 
there  exist  two  interrelated  aspects:  (i)  Find  a  method  for  generating 
particular  solutions.  ( i i )  Find  a  method  for  solving  the  initial  value 
problem.  With  respect  to  (i)  we  note:  (a)  There  exist  certain  explicit 
transformations  which  map  solutions  of  a  given  Painleve  equation  to  solu¬ 
tions  of  the  same  equation  but  with  different  values  of  the  parameters. 
Such  maps,  called  Schlesinger  transformations ,  were  given  by  Lukashevich 
and  Gromak  [31]-[34]  for  PII-PV  and  by  Fokas  and  Yortsos  [35]  for  PVI. 
Furthermore  it  turns  out  that  for  certain  choices  of  their  parameters  PII- 
PV  I  admit  rational  solutions  as  well  as  one-parameter  families  of  solu¬ 
tions  expressible  in  terms  of  Airy  [2],  [36],  Bessel  [37],  Weber-Hermi te 
[38],  Whittaker  [39]  and  hypergeometri c  [40]  functions  respectively. 

Using  the  above  transformations  and  special  solutions,  one  can  construct 
(for  certain  choices  of  the  parameters)  various  elementary  solutions  of 
PII-PVI.  These  solutions  are  either  rational  or  are  functions  which  are 
related,  through  repeated  di fferentiations  and  multiplications,  to  the 
above  mentioned  classical  transcendental  functions,  (s)  Ablowitz  and 
Segur  [10]  characterized  a  non-elementary  one  parameter  family  of  solu¬ 
tions  of  P 1 1  through  a  Gel 1 fand-Levi tan-Marchenko  integral  equation  of 
the  Fredholm  type,  (y)  Fokas  and  Ablowitz  [41]  characterized  a  two- 
parameter  family  of  solutions  of  P 1 1  using  a  matrix  system  of  Fredholm 
integral  equations.  However,  in  both  (e)  and  (y)  the  free  parameters  were 
not  related  to  the  initial  data  of  PII. 


The  main  focus  of  this  paper  is  to  give  a  method  to  solve  the 
initial  value  problem  of  the  Painleve  equations.  In  this  respect  we 
note:  (i)  Flaschka,  Newell,  Jimbo,  Miwa  and  Ueno  introduced  a  new 
powerful  approach  for  studying  the  associated  initial  value  problem; 
solving  such  an  initial  value  problem  is  essentially  equivalent  to  solving  an 
inverse  problem  for  a  certain  isomonodromic  linear  equation  (see  (iv) 
above).  In  analogy  with  the  1ST  method  introduced  by  Gardner,  Greene, 

Kruskal  and  Miura  [42]  we  call  the  above  method  an  inverse  monodromic 
transform  method  (IMT).  (ii)  Flaschka  and  Newell  [14]  applied  the  above 
method  to  the  solution  of  PI  I  and  to  a  special  case  of  Pill.  They  form¬ 
ulated  the  inverse  problem  in  terms  of  what  the  authors 
of  [14]  call  a  system  of  singular  integral  equations,  (iii)  Jimbo,  Miwa 
and  Ueno  [15]  considered  the  Painleve  equations  within  the  larger  program 
of  study  of  monodromy  preserving  deformations  for  a  first  order  matrix 
system  of  ODE's  having  regular  or  irregular  singularities  of  arbitrary 
rank.  The  inverse  problem  is  solved  in  terms  of  formal  infinite  series 
uniquely  determined  in  terms  of  certain  monodromy  data,  (iv)  Fokas  and 
Ablowitz  demonstrated  that  the  inverse  problem  of  PI  I  can  be  formulated 
as  a  matrix,  singular,  discontinuous,  homogeneous  Riemann-Hi 1 bert  (RH) 
problem  defined  on  a  complicated  contour.  This  has  conceptual  and  prac¬ 
tical  implications:  Conceptually,  it  becomes  clear  that  there  is  a 
unified  approach  to  solving  certain  initial  value  problems  for  equations  in  1, 
1+1  (one  spatial  and  one  temporal)  and  2+1  dimensions.  Using 
techniques  from  RH  theory,  the  RH  problem  can  be  simplified  sub¬ 
stantially  (it  can  be  mapped  to  a  series  of  regular,  continuous  RH  problems, 
each  defined  on  the  real  axis). 
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In  this  paper  we  present  a  general  method  for  solving  the  initial 
value  problem  associated  with  a  given  Painleve  equation.  This  method, 
which  simplifies  and  extends  ideas  of  [43],  involves  three  main  steps: 

1.  Use  classical  theory  of  linear  ODE’s  to  formulate  a  RH  problem  for  a  func¬ 
tion  called  f(z,t)  which  solves  the  underlying  isomonodromic  linear  equa¬ 
tion.  This  basic  RH  problem  is, in  general,  a  matrix  singular,  discontin¬ 
uous,  problem  formulated  on  a  complicated  contour  (several  intersecting 
rays)  and  is  uniquely  defined  in  terms  of  certain  monodromy  data.  2. 

Choose  the  parameters  of  the  given  Painleve  equations  in  such  a  way  that  the  above 

problem  is  nonsingular;  then  map  the  basic  RH  problem  to  series  of  RH  prob¬ 

lems  defined  on  simple  contours.  All  of  these  problems  except  one  are 
continuous.  Furthermore,  some  of  them  can  be  solved  in  closed  form  (in 
terms  of  a  quadrature).  Use  certain  auxiliary  functions  to  map  the  dis¬ 
continuous  RH  problem  to  a  continuous  one.  Then  apply  the  rigorous  re¬ 
sults  of  the  RH  theory,  e.g.  [44]  to  establish  the  existence  and  unique¬ 
ness  of  the  solutions  of  the  above  continuous  RH  problems.  3.  Use  the 
basic  RH  problem  to  obtain  Schlesinger  [51],  [13],  transformations,  shifting  by 
an  integer  or  by  a  half  integer  all  the  parameters  of  the  given  Painleve 

equation.  Hence,  using  these  transformations,  the  study  of  the 

singular  RH  problem  reduces  to  the  study  of  the  regular  one.  We  note 
that  for  special  choices  of  the  monodromy  data  the  basic  problem  can  be 
solved  in  closed  form.  This  yields  particular  solutions  of  Painleve 
equations  expressible  in  terms  of  the  classical  transcendental  functions 
mentioned  above. 

The  above  method  is  applied  to  the  solution  of  the  initial  value 
problem  of  PIV  and  PV:  The  RH  problem  corresponding  to  each  of  these 
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equations  is  mapped  to  two  RH  problems  on  simple  contours,  one  of  which 
can  be  solved  in  closed  form,  while  the  other  can  be  made  continuous. 

Furthermore,  Schlesinger  transformations  are  derived  for  both  PIV  and 
PV.  For  special  choices  of  the  monodromy  data  the  basic  RH  problems  for 
PIV  and  PV  can  be  solved  in  closed  form;  this  yields  solutions  of  PIV 
and  PV  in  terms  of  Weber-Hermi te  and  Whittaker  functions  respectively. 

PII  is  considered  in  [45],  using  a  different  i somonodromic  spectral 
problem  than  the  one  used  in  [43].  These  results,  and  Pill  (which  is 
related  to  a  special  case  of  PV  [9])  will  be  presented  elsewhere.  PVI 
has  been  solved  by  C.  Cosgrove  [13]  and  PI  remains  open. 

The  Hamiltonian  structure  of  the  Painleve  equation  is  studied  in 
[14]  and  [56]. 

2.  THE  GENERAL  FRAMEWORK 
2.1.  RH  Problems. 

Let  C  be  a  simple,  smooth,  closed  (or  infinite)  contour  dividing  the  complex  z- 
plane  into  two  regions  D+  and  D~  (the  positive  direction  of  C  will  be 
taken  as  that  for  which  D+  is  on  the  left). 

0* 


Figure  2.1 


A  function  $(z\  defined  in  the  entire  plane,  except  for  points  on 
C  which  will  be  called  sectionally  holomorphic  if:  i)  the  function  $(z) 


is  holomorphic  in  each  of  regions  D+  and  D~  except,  perhaps,  at  2  =  °° ; 
ii)  the  function  #(z)  is  sectionally  continuous  with  respect  to  C, 
approaching  the  definite  limiting  values  $+(c),  $'(0  as  2  approaches 
a  point  c  on  C  from  D+,  or  O',  respectively.  The  classical  homogeneous 
RH  problem  is  defined  as  follows  [46].  Given  a  contour  C,  and  a  function 
G(c)  which  is  Holder  on  C  and  det  G(t)  t  0  on  C,  find  a  sectionally 
holomorphic  function  i>(z),  with  finite  degree  at  such  that 

$+U)  =  GU)$"U),  on  C,  (2.1) 

where  $  (;)  are  the  boundary  values  of  t>(z)  on  C.  If  G(c)  is  scalar,  (2.1) 

able  in  terms  of  quadratures.  If  G(c)  is  a  matrix  valued  function, 

then  (2.1)  is  in  general  solvable  in  terms  of  a  system  of  Fredholm 

integral  equations.  Various  generalizations  of  the  above  RH  problem  are 

possible.  For  example:  i)  The  contour  C  may  be  replaced  by  a  union  of 
* 

intersecting  contours,  ii)  G(c)  may  have  simple  discontinuities  at  a 
finite  number  of  points;  in  this  case  one  allows  i(z)  to  have  mtegrable 
singularities  in  the  neighborhood  of  these  points,  iii)  RH  problems  may 
be  considered  in  other  than  Holder  spaces  (e.g.[47]):  iv)  One  may  con¬ 
sider  inhomogeneous  RH  problems  ♦*(;)  =  G U)$*(c)  +  F(c)  on  C. 

It  is  interesting  that  the  first  RH  problem  was  formulated  in  con¬ 
nection  with  an  inverse  problem  (see  [43]  for  references).  Actually,  RH 
problems  are  intimately  related  to  the  solution  of  inverse  problems  in 
1+1  (one  spatial  and  one  temporal),  2+1  and  1  dimensions: 

2.2.  Inverse  Problems  in  1+1 

We  recall  that  a  necessary  condition  for  a  given  nonlinear  equation 
for  q(x,t)  to  be  solvable  via  1ST  is  that  this  equation  is  the  compati- 


is  so 


bility  condition  of  a  pair  of  linear  equations.  Let  us  consider  the 
modified  KdV 


^  *  Vx  -  * 0  <2-2' 

as  an  illustrative  example  [48].  Equation  (2.2)  is  the  compatibility 
condition  of 


i*x(c.x,t)  =  c 


H'(t,X,t)  + 


q(x,t) 


,q(x,t)  0 


f(;,x,t). 


(2.3a) 


^t(  C  »x ,  t) 


/  -4ic3-2iq2c 
\4q42-2iqxt-2q3-qxx 


4qc2+2iqxc+2q3-q 

3  2 

4i;J+2iq  c 


f(c,x,t). (2.3b) 


We  first  note  that  the  above  Lax  pair  is  isospectral ,  i.e.  dc/dt  =  0. 

Also  it  turns  out  that  equation  (2.3a)  is  of  primary  importance;  equation 
(2.3b)  plays  only  an  auxiliary  role.  To  solve  the  initial  value  problem, 
for  initial  data  decaying  as  |x|  -  »  ,  one  first  formulates  an  inverse 
problem  for  f(z,x,t):  Given  appropriate  scattering  data  reconstruct  't . 

By  studying  the  analytic  properties  of  f  with  respect  to  z,  where  f  sat¬ 
isfies  (2.3a),  one  establishes  that  there  exists  a  f  which  is  a  section- 
ally  meromorphic  function  of  z,  with  a  jump  along  the  Re  z  axis.  This 
jump,  as  well  as  the  residues  of  the  poles,  are  given  in  terms  of  approp¬ 
riate  scattering  data.  Thus  the  inverse  problem  is  equivalent  to  a 
matrix,  regular,  continuous,  RH  problem  defined  along  the  Re  z  axis  and 
uniquely  specified  in  terms  of  the  scattering  data. 

Since  in  the  above  discussion  we  have  only  used  (2.3a),  it  is 


evident  that  one  may  pose  an  inverse  problem  for  an  appropriate  function  q(x).  How- 


-9- 


ever,  this  result  is  useful  for  solving  the  initial  value  problem  of 
q(x,t)  only  if  q  evolves  in  such  a  way  in  t,  that  the  scattering  data  is 
known  for  all  t.  If  4»  evolves  in  t  according  to  (2.3b)  (i.e.  if  q 
solves  (2.2))  then  it  turns  out  that  the  evolution  of  the  scattering 
data  with  respect  to  t  is  simple.  Hence  the  above  RH  problem  is  speci¬ 
fied  in  terms  of  initial  scattering  data;  its  solution  yields  v(c,x,t) 
and  then  (2.3a)  gives  q(x,t). 


2.3.  Inverse  Problems  in  2+1. 

Let  us  consider  the  Davey-Stewartson  equation  (a  two  dimensional 
analogue  of  the  nonlinear  Schrodinger  equation) 


*  ?Kx*V  =  -X|0|2|3n0-  ♦xxAy'2*(IQ|Z>XXi  '  =  ° =  =1 


(2.4) 


as  an  illustrative  example  [49].  A  Lax  pair  for  (2.4)  is  given  by 

10  f 0  °\ 

Y  =  i;(J»-yJ)  +  qt  +  oJt  ,  Jr  (i  ),  q?L  )! 

*  y  u  -i  \^xQ  0/ 


(2.5a) 


't  =  Vyy  +  A2fy  +  Alf  "  ( A3'*'“yA3o)  +  2UA3fy  +  (2.5b) 


where  A^,  A^ ,  A^,  A^q  are  appropriate  matrix  functions  of  Q,  Q  (  Q  denotes 
the  complex  conjugate  of  Q). 

The  situation  is  conceptually  similar  to  the  case  of  1+1:  To  solve 
the  initial  vlaue  problem  for  q(x,y,t)  one  first  formulates  an  inverse 
problem  for  f(z,x,y,t).  Depending  on  the  value  of  a  there  exist  two 


different  cases  (for  brevity  of  presentation  we  assume  non-existence  of 
poles,  i.e.  non-existence  of  lumps):  (i)  a  =  -1.  There  exists  a  ?  which 
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is  a  sectionally  holomorphic  function  of  z  and  which  has  a  jump  along 
the  Re  z  axis.  This  jump  is  also  given  in  terms  of  scattering  data  but 
it  depends  on  them  in  a  non-local  way.  Thus  the  inverse  problem  is 
equivalent  to  a  non-local ,  matrix,  regular,  continuous  RH  problem  defined 
along  the  Re  z  axis  and  uniquely  specified  in  terms  of  scattering  data. 

(ii)  o  «  i.  There  exists  a  f  which  is  bounded  for  all  complex  z,  but 
which  is  anlaytic  nowhere  in  the  complex  z  plane.  However,  its  departure  from 
holomorphici ty  3H73z  can  be  expressed  linearly  in  terms  of  f  and  appropriate 
inverse  data.  Thus,  now  the  inverse  problem  is  equivalent  to  a  3  (DBAR) 
problem:  Given  3V/3Z  reconstruct  f.  The  3  problem  is  a  generalization  of 
a  RH  problem  and  has  been  studied  extensively  in  the  mathematical  lit¬ 
erature  [50]. 

Using  (2.5b),  again  one  shows  that  the  inverse  scattering  and  the 
inverse  data  evolve  simply  in  time.  Hence,  the  above  RH  and  3  problems 
are  specified  in  terms  of  initial  data;  their  solutions  yield  1*'(c,x,t,t) 
and  then  (2.5a)  gives  q(x,y,t). 


2.4,  Inverse  Problems  in  0+1. 

The  Lax  pair  associated  with  the  PIV  equation 


d  y  .  l_  ,dv,2 


♦I 


4ty2  ♦  2(t2  ♦  a)y  ♦  £ 


is  given  by  [15] 


(2.6) 


Y 

Z 


(z) 


j .  ■(', 

•v  Vu(u-ve-> 


V' 

— (v-2e  ’ 
.uy  0' 


HI 

2 


(2.7a) 


1 
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v ( 2 ) .  (2.7b) 


V2>  = 


2  + 


0  -1 


V'-’o-V 


Indeed  Y 


xt 


Y tx  implies 


at  =  -4*  ♦  y2  *  2ty  *  %•  af  *  -u<y*2t)’ 


40, 


I  =  -  f  V2  +  M  -  yj  v  ♦  (e0  +  Qjy, 


(2.8) 


where , 


a  =  20  -  1,  0  =  -86* 


(2.9) 


As  in  the  cases  of  1+1  and  2+1,  solving  the  initial  value  problem 
of  PIV  reduces  to  solving  an  inverse  problem  for  Y:  Reconstruct  Y(z,t) 
in  terms  of  appropriate  monodromy  data.  Again  this  inverse  problem 
will  be  solved  in  terms  of  a  RH  problem.  Thus  it  is  essential  to  study 
the  analytic  properties  of  Y  with  respect  to  z.  However,  in  contrast  to 
the  analogous  1ST  problem  in  1+1  and  2+1,  the  task  here  is  straight¬ 
forward:  Equation  (2.7a)  is  a  linear  ODE  in  z,  therefore  its  analytic 
structure  is  completely  determined  by  its  singular  points.  In  this  part¬ 
icular  case  z  =  0  is  a  regular  singular  point  and  z  =  °»  is  an  irregular 
singular  point  of  rank  2.  Complete  information  about  the  singular  point 
z  =  0  is  provided  by  the  monodromy  matrix  Mg.  Complete  information  about 
z  =  «  is  provided  by  the  monodromy  matrix  and  by  the  Stokes  multipli¬ 
ers  a,  b,  c,  d.  Solutions  of  (2.7a),  Yg  and  Y ^ ,  normalized  at  zero  and 
infinity  respectively  are  related  via  a  connection  matrix  Eg  with  entries 
ag,  Bg,  Yg,  <5g.  Taking  into  consideration  the  above  singularities,  there 
exist  a  sectionally  holomorphic  function  Y,  with  jumps  across  four 
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rays,  arg  z  s  -  |,  and  with  singularities  at  z  =  0,  z  =  *.  The 

jumps  are  specified  by  a,  b,  c,  d  and  the  nature  of  singularities  by  MQ, 
M^.  This  leads  to  a  matrix,  singular,  discontinuous  RH  problem,  defined 
on  the  above  rays  and  specified  in  terms  of  the  monodromy  data 

Monodromy  Data  (MD)  =  {a,  b,  c,  d,  aQ,  SQ,  yq,  6q; .  (2.10) 

Consistency  of  the  above  RH  problem  yields 


4 

(  n 
j  =  l 


Gj)M-  = 


(2.11) 


where  Gj  are  the  Stoke  matrices  uniquely  defined  in  terms  of  the  Stokes 
multipliers.  Using  (2.11)  and  certain  similarity  arguments  it  can  be 
shown  that  all  MD  can  be  expressed  in  terms  of  two  of  them.  Furthermore, 
x  equation  (2.7b)  implies  that  the  MD  are  time-invariant.  Hence  the  above 
basic  RH  is  specified  in  terms  of  two  initial  parameters  (these  two  initial 
parameters  are  obtained  from  the  two  initial  data  of  PIV).  The  solution 
of  this  RH  problem  yields  Y(z,t)  and  hence  (2.7a)  yields  y(t). 

This  RH  problem  can  be  simplified  considerably :  (i) 

Assume  0  <  <  1 ,  0  1  e0D  <  1 »  9q  t  then  the  above  RH  problem  is  non¬ 

singular.  It  is  interesting  that  the  basic  RH  problem  can  be  used  to  obtain 
Schlesinger  transformations  which  shift  6g  and  by  a  half-integer. 

By  using  these  transformations  the  general  case  is  reduced  to  the  regular 
case,  (ii)  The  basic  RH  problem  can  be  mapped  to  a  sequence  of  two  RH 
problems,  one  on  the  line  arg  z  =  4  and  the  other  on  the  line  arg  z  =  -  ^. 
The  first  one  is  continuous  (both  at  z  =  0  and  z  =  ®);  furthermore  it 


can  be  solved  in  closed  form.  The  second  one  is  discontinuous  both  at 
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2=0  and  2  =  ».  By  using  standard  auxiliary  functions  [46]  this  discon¬ 
tinuous  problem  is  mapped  to  a  continuous  one.  The  theory  of  con¬ 
tinuous  RH  problems  on  simple  contours  can  then  be  used  to  establish  uniqueness 
and  existence  of  solutions.  Elementary  solutions  of  PIV,  expressible  in 
terms  of  Weber-Hermi te  functions  are  obtained  naturally  within  the  above 
formalism. 

We  hope  that  the  above  discussion  elucidates  the  connection  between 
1ST  (inverse  scattering  transform)  and  IMT  (inverse  monodromic  transform): 

There  exists  a  unified  approach  to  initial  value  problems  in  1,  l+l,  and 
2+1  dimensions:  Solving  the  initial  value  problem  of  an  integrable  equation: 
q(t)  or  q(x,t)  or  q(x,y,t)  is  equivalent  to  solving  an  inverse  problem 
fora  suitable  eigenfunction  f(z;t)  or  Y(z;x,t)  or  y(z;x,y,t).  The 
inverse  problem  generically  takes  the  form  of  a  RH  problem  for  equations 
in  1,  1+1,  and  in  general  the  form  of  a  3  ( DBAR)  problem  for  equations  in 
2+1  (the  DBAR  problem  being  a  general  Nation  of  a  RH  problem).  To  define  the 
relevant  RH  or  DBAR  problems  one  needs  to  study  the  analyticity  properties 
of  f  with  respect  to  z.  Furthermore  these  problems  are  uniquely  defined 
in  terms  of  certain  asymptotic  data  of  the  underlying  linear  system  satis¬ 
fied  by  Y  (monodromy  data  in  the  case  of  equations  in  1  dimension  and 
scattering  data  in  the  case  of  equations  in  1+1  and  2+1). 

We  note  that  the  linear  limit  of  the  1ST  yields  the  Fourier  trans¬ 
form  of  q(x,t).  In  that  sense  1ST  is  the  nonlinear  analogue  of  the  Fourier 
transform  [52].  Since  the  linear  limit  of  the  IMT  is  the  Laplace's  method 
for  linear  ODE's,  the  IMT  is  the  nonlinear  analogue  of  the  Laplace's 


method. 
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2.5.  Solution  of  a  Matrix  Continuous  RH  Problem. 

From  the  above  discussion  it  follows  that  solving  the  initial  value 
problem  of  a  Painleve  equation  reduces  to  solving  a  matrix  continuous  RH 
problem  along  a  simple  contour.  We  recall  that  the  solution  of  the  RH  problem 
(2.1),  where  G(c)  is  Holder  (i.e.  all  its  entries  satisfy 
|  G  j  k  ( C  i )  -  G  j  k  ( C  2 )  I  <  A I  ^  i  2  ^  A  *  for  some  constants  A  an<*  0  <  A  1  1,  for  all 
?  on  C)  and  det  G(c)  t  0  on  C  is  given  by 


♦'«>  -  m 


dc[G~1(4)G(;)-I]<t(c)  = 

C  c  -  c 


(2.12) 


where  ♦*  is  the  value  of  <t(z)  at  infinity  (we  assume  that  <t  has  a  finite 
degree  at  infinity).  Equation  (2.12)  can  be  obtained  by  writing  the 
conditions  that  <t+,  are  +  and  -  functions  respectively  and  then  replac 
ing  <t>+  by  G<t~. 

In  what  follows  we  define  RH  problems  on  suitable  rays.  These  rays 
are  naturally  defined  for  a  given  problem  (see  §3.2). 


i 
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3.  PAINLEVE  IV 

In  this  section  we  consider  the  fourth  Painleve  equation  (2.6). 

We  first  use  equations  (2.7)  to  study  the  analytic  properties  of  Y(z,  t), 
as  well  as  the  properties  of  the  monodromy  data. 


3.1.  The  Direct  Problem. 

Proposition  3.1. 

Let  Yq  be  the  solution  of  (2.7a)  analytic  in  the  neighborhood  of 

2=0  and  normalized  by  the  requirements  that  det  Yq  =  1  and  that  YQ  also 

solves  (2.7b).  Let  Y.,  j  =  1,...,4  be  solutions  of  (2.7a)  analytic  in 

the  neighborhood  of  infinity  such  that  det  Y^  =  1  and  Y  ^  as  jx|  -  °° 

in  S-,  Y  is  the  formal  solution  matrix  of 
J 

(2.7a)  in  the  neighborhood  of  infinity,  and  the  sectors  Sj  are  given  by 


Figure  3.1 


Then  the  analytic  functions  Yg,  Y^,...,Y^  satisfy: 

(i)  Yq( z )  x  YQ( 2 ) 2  0  as  z  -  0;  Dq  *  Diag(e0,-e0) ,  6Q  f  n 


where  Yg(z)  is  holomorphic  at  z  =  0.  (If  0g  =  n/2,  Yg(z)  has 
a  logarithmic  singularity). 

(ii)  Y  .(z)  %  Yjzje^2^  ^  )  as  \z\  -  •,  z  in  S- ,  *  Diag(e 


Q(z)  t  Diag(q.-q),  q{z,t)  *  f-  +  zt,  Y  (z)  is  holomorphic  at 


z  =  «*,  and  Y  (z)  I  +  0(4) . 


(iii)  YQ(zez1 " )  =  YQ{ z )MQ ,  MQ  * 


2  i  TT0n  2inen 

e  u  2inJe  u 

-2irre_ 


,  (3 


J  =  0  if  e0  ,  3*1  if  eQ 


(iv)  YJz)  =  Y.  (  z  )G,  ,  Y,(z)  =  Y„(  z  )G0  ,  Y.(Z)  »  Y,(z)G,, 


w 


Furthermore,  the  parameters 


MD  *  {a,b,c,d,a0,8g,y0,60} 


satisfy  the  following  consistency  condition. 


(vi) 


V"-  ■  ESVEo- 


(3.8) 


(3.9) 


Proof. 

1.  Analysis  near  2=0: 

It  is  well  known (see  for  example  [53]) that  if  the  coefficient  matrix  of  a  linea 
differential  equation  has  an  isolated  singularity  at  z  -  0,  the  solution  of 
the  differential  equation  will  in  general  he  singular  at  z  =  0. 

This  solution  can  be  obtained  in  the  form  of  a  formal  power  series;  this 
series  actually  is  convergent  in  an  appropriate  circle  of  the  com¬ 
plex  plane.  In  this  particular  case  if  YQ  =  (Y^,  Y^)  we  find: 


(3.10a) 


(z) 


where , 


H(|)  =  -y(v-0Q-eJ  -  vt  *  (1  +  29n-v)(t  ♦  2vy), 
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"Si!  ♦  f  [t  ♦  ^  -  rzvTi  <v  -  %  ■ 


c(t)  *  jV 


(Expressions  for  H ^ 2 j may  a^s0  given,  but  are  not  necessary  for  our 
discussion).  The  multiplicative  constants  with  respect  to  z  in  (3.10)  are  fixed 
by  the  requirement  that  (3.10)  also  satisfy  (2.7b).  We  note  that  when 
Oq  =  n/2  there  will  be,  in  general,  a  logarithmic  term  and  the  two 
linearly  independent  solutions  have  the  form: 


Yq^(z)  ,  Yq2)(z)  =  J(zn  z)Y^1)  +  z  0  Y^2)(z),  Oq  =  \  ,  (3.11) 

-  (2) 

where  is  a  polynomial  in  z  and  J  is  a  complex  constant.  Equations 

21 ’Dn 

(3.10)  imply  Y g(ze^  n)  =  Yg(z)e  .  Similarly,  equations  (3.10)  and 

(3.11)  imply  (3.4). 


2.  Analysis  near  z  =  «. 

The  two  linearly  independent  formal  solutions  Y  (z)  - 

1 2 1 

Y '  ' ( z ) )  of  (2.7a)  have  the  expansions: 


( fl1 ’ ( z) , 


xiU(z)  = 


y<2)(z)  . 


;i)  V 


(3.12a) 


(3.12b; 


where 


K  !  J(,  -  20o)  -  (v  -  0O  -  0j(t  *  X). 

Suppose  that  Y.  is  the  asymptotic  expansion  of  Y^  for  large  z  in  S, 


According 


to  the  Stokes  phenomenon,  the  asymptotic  expansion  of  Y^  in  sector  S£ 
is  given  by  Y  6  where  G  is  a  constant  matrix.  A1 ternati vely,  one  may 


introduce  different  solutions  Yj,...,Y4  such  that  Y  ^  is  asymptotic  to 

Y  in  S..  Then,  for  example,  since  both  Y. ,  Y,  solve  (2.7a)  it  follows 
00  j  i  <- 

that  Y ^  =  Y^Gp  where  G^  is  a  constant  (with  respect  to  z )  non-singular  matrix;  using 
the  asymptotic  expansions  of  Y^,  Y ^  it  can  be  shown  that  G^  has  the  form  given  in  (3.6 
[54].  Similarly  we  take  Yg'v  Y^  in  S^,  where  S^:  2*  -  ^  <  arg  z  <  2,  + 
then  Y,.  *  Y4G4.  Birkof  [54]  has  related  Y^  and  Y ^ :  Y^  and  Y,.  are  defined 
for  all  z,  however,  they  both  tend  to  Y^  only  in  S,  and  S^  respecti vely: 


YjUj)  -vYJzjJe 


«ll> A 
Z1 


VZ5' 


0(2*;)  ,  1 
le  b  (f- ) 
Z5 


If  z^  is  in  Sj  then  z^e2’"  is  in  S^.  Thus 


0(2,  )  ,  0  -2i^D 


Q(z,e2’’T)  .  D  Q(z.)  ,  D.  ■ 

Ys(2ie  >  '-Y-<2ie  >e  (-Tn)  =  Y.<2 i>*  <rr)  e 

Z  le  1 


where  we  have  used  that  Y^z)  is  holomorphic  at  z=°°  and  Q(z^e^  ")  r  Q(z^).  Thus 


Y 1  ( z )  -  Y 5(ze2l,T)e 


(3.13) 


Hence 

o  •  -•  »■  23^0 

Y 5(z)  =  Y4(z)G4,  or  Y5(ze2  )  «  Y4(ze^  )G4,  or  Y^z)  =  \{ze'  )G4e 


3.  Connecting  Yq  and  Y ^ . 

Since  both  YQ  and  Y^  satisfy  (2.7a),  Y^  =  YqEq  and  det  EQ  -  1, 
since  det  Yj  =  det  YQ  = 


1.  Eq  is  called  a  connection  matrix. 
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4.  Consistency. 

We  first  note  that  if  the  solution  matrices  ,  Y^  are  related  via 
the  constant  matrix  C,  and  if  M^,  are  their  monodromy  matrices  about 
the  same  regular  singular  point  then 

Y1  =  Y2C  —  ”  M1  =  c'1m2C  ;  (3.14) 


this  because: 


2i  n 


2i  TT 


YjUe  )  =  Y,(  ze  )C  «  Y,(z)M,C  «  Y,C  %C. 


2 v * 1  2 


Equations  (3.5)  imply  Y ^ ( z )  =  Y4(ze2i ’Mg^  =  Y3(ze2i n)G3G4Ma 


=  Y  ,(ze411t)(  n  G  ■ )  M_ .  But  since  Mn  is  the  monodromy  matrix  of  YQ  and 


j  =  l 


J 


■-U 


Yi  =  YoEo*  then  Eo  MqEq  is  the  monodromy  matrix  of  Yj.  Thus 

,  4 

Y ^ ( z )  =  Y^(z)Eq1MqEq(  n  GjM^,  which  implies  (3.9). 

jsl 


Remark  3.1. 

(i)  One  has  two  choices:  either  to  consider  four  different  fundamental 

solutions  Y Y.  such  that  Y.  %  Y  as  f  Z i  -  =°,  z  in  S-,  or  to  con- 

sider  one  fundamental  solution  Y,  but  then  Y  a.  G.Y  as  i  z 1  *  »  ,  z  in  S,. 

J  00  J 

i 

We  intend  to  use  these  Y's  to  formulate  a  RH  problem,  hence  it  is  lmport- 

J 

ant  to  have  the  same  behavior  at  infinity,  that  is  why  we  choose  four 
different  solutions. 

I 

(ii)  The  solutions  Y^,...,Y4  are  defined  in  the  whole  complex  z  plane  and 
the  relationships  (3.5)  are  valid  everywhere.  However,  in  order  to  form¬ 
ulate  a  RH  problem  we  will  restrict  the  domain  of  the  Y^'s  only  in  the  sector  t 
that  their  asymptotic  behavior  is  Y^,  thus  we  will  use  (3.5)  only  on  the 


! 
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rays  C1,...,C4  (see  Figure  3.1). 

Proposition  3.2. 

(i)  The  monodromy  data,  MO,  given  by  (3.8)  and  defined  in  Proposition 
3.1  are  time-invariant. 


(ii)  All  of  the  MO  can  be  expressed  in  terms  of  two  of  them.  This 
follows  from: 

1.  det  Eq  =  1. 

2.  Equation  (3.9). 

3.  If  Y  solves  (2.7)  with  y  satisfying  PIV,  then  Y  r  R  *YR, 

R  r  diag(r1^,  r"1^),  where  r  is  an  arbitrary  complex  constant, 
also  solves  (2.7)  with  y  satisfying  PIV.  The  Stokes  matrices  Gj 
and  the  connection  matrix  Eq  are  transformed  to  Gj  *  R  *GjR, 

E„  !  R-^R,  i.e. 

a  =  ra,  b  =  b/r,  c  =  rc,  d  =  d/r,  =  Xq,  =  eQ/r, 


r0 


rY  o  *  =  80* 


(3.15) 


Thus  r  may  be  chosen  to  eliminate  one  parameter,  e.g.  r  =  8Q. 

4.  Changing  the  arbitrary  integration  constant  of  o(t)  (see 


(3.10)  )  amounts  to  multiplying  Y^(z),  Y^(z)  by  the  arbitrary 


,(2) 


complex  constants  p  and  p’^  respectively.  This  maps  Eq  t< 


L  *  PEn,  P  *  Diagjp.p'1) ,  i  .e. 


PaQ ’  80  =  PV  Yo 


y_0  :  _ 0 
p  *  0  '  P 


(3.16! 


Thus  p  may  be  chosen  to  eliminate  one  parameter,  e.g.  p  -  .q. 


lii)  Equation  (3.9)  implies 
2  i  rrG 

(l+bc)e  “  *  [ad  ♦  ( 1+cd) ( l+ab) ]e 


2cos2-!-  .  (3.17 
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Proof . 

(i)  Similar  to  the  proof  given  in  [14]. 

(ii)  Y  solves  (2.7)  iff  y  =  y,  v  =  v,  u  =  ur"^  which  is  consistent  with 

(2.8).  Parts  (ii)  3  and  (iii)  are  straightforward.  Equations  (3.15), 

(3.16)  may  be  chosen  to  fix  two  of  the  entries  of  the 

connection  matrix.  Then  det  Eg  =  1  and  equation  (3.9)  imply  the  rest 
of  the  MO  in  terms  of  two  of  them. 

3.2.  The  Inverse  Problem. 

In  what  follows  we  formulate  a  RH  problem  for  the  case  that 

0  <.  0g  <  1 ,  0  <  0^  <  1.  This  assumption  leads  to  a  regular  RH  problem. 

The  general  case  follows  by  considering  the  results  of  this  section  and 
of  §3.3. 

In  what  follows  we  shall,  for  convenience  of  notation,  consider  RH 
problems  along  suitable  rays.  Actually,  in  order  to  satisfy  convergence 
criteria  these  rays  must  be  deformed  appropriately.  The  deformation  process 
is  to  connect  at  large  values  of  z  the  rays  to  "asymptotic"  curves  defined 
by  Re  q(z,t)  *  0.  These  asymptotic  curves  tend  to  the  straight  line  rays 
(i.e.  Re  q(z,t=0))  for  |zj  ♦  ®. 

An  alternative  procedure  which  we  anticipate  to  be  equivalent  (but 
one  which  we  have  not  seriously  considered)  is  to  deform  the  rays  by  a 
sufficiently  small  angle  e  into  the  region  Re  q(z,t)  <  0.  In  this  case  the 
RH  problem  has  jump  matrices  which  rapidly  tend  to  unity  as  z  *  ».  We 
expect  the  limit  as  €  +  0  of  this  deformed  RH  problem  should  tend  to  the 
solution  of  the  RH  problem  discussed  above. 

Theorem  3.1. 

Consider  the  following  matrix,  regular,  homogeneous  RH  problem  along 

the  four  rays  Cp...,C^  (Figure  3.1):  Determine  the  sectionally  holo- 

morphic  function  t(z),  *(z)  =  vAz)  if  z  is  in  S  ,  j  =  1,...,4,  from  the 

J  J 

following  conditions: 


1. 


v  ■  satisfy  the  jump  conditions 

O 


f2(c)  =  t1(c)g1(c),  f3(c)  =  f2(c)g2(c),  »4U)  =  f3(;)g3U), 

=  4'4(ce2l,T)g4(0  (3.18) 

along  the  rays  C2,  C3>  C4,  respectively,  where 

9J  #  eQGj6"Q’  J  =  1,2’3>  g4  *  eQG4e_QM»-  (3.19) 

2.  ^(z)  %  (^)  (I  +  0(-))  as  |z|  00  *  in  S..  (3.20) 

3.  f(z)  has  at  most  an  integrate  singularity  at  the  origin  with  a 
monodromy  matrix  given  by 

»1(ze21ir)  =  ^(zJE^MqEq,  z  -  0.  (3.21) 


In  the  above,  G j ,  Q,  M^,  0^,  Mq  are  defined  in  Proposition  3.1. 


4.  The  monodromy  data  MD,  given  by  (3.8),  satisfy  the  properties 
given  in  Proposition  3.2  (ii).  Then: 

(i)  The  above  RH  problem  is  discontinuous  both  at  the  origin  and  at 
infinity.  Actually 


(3.22) 

f o 1  low¬ 
ing  RH  problem  along  the  contour  C1  ♦  C3:  Determine  the 


0;  n  g  .  \  M  ,  z  -  *  . 
j  =  l  J 

(ii)  To  obtain  the  solution  of  the  above  RH  problem  consider  the 


Figure  3.2 


secti ona  1 1  y  holomorphic  function  k(z),  k(z)  =  k,  (z)  if  z  is  in  S  +  <; 

i  J  ^  J2  * 

k(z)  =  k2(z)  if  z  is  in  +  S^,  from  the  following  conditions: 


k.  satisfy  the  jump  condition 


k  *  k 
*1  *2 


■1  -be 


y  ■ 


on  Cj, 


h(z)  * 


ao(z)  1 


°<z>  *  ‘  ST 


,-2qU) 


(If  h^,  h^  denote  h  in  and  S2  +  respectively  then  h 

h  =  h2  on  C3) . 

,  0  . 

2.  k(z)  %  (|)  *( I  ♦  0(1))  as  | z 1  -  •  . 


(3.23) 


—  h^  on  C| j 


(3.24) 


3.  k(z)  has  at  most  an  integrable  singularity  at  the  origin  with  a 


monodromy  matrix  given  by 


k(ze2^)  =  k(z)h1(0)EQ1MQE0h'11(0),  z  -  0. 


(3.25) 


The  above  RH  problem  is  discontinuous  both  at  the  origin  and  at  infinity. 

Actually  if  g.  g.  denote  the  jump  matrices  along  C, ,  C,  respectively 
*  1  j  1  J 


-  ” 1 u“ 1  c  u-1 


%l\-  hl<°’E0  M0  EOhl 


,  Z  *  x. 


(3.26) 


However,  the  above  RH  problem  can  be  mapped  to  a  continuous  one  using  the 


-25- 


the  auxiliary  functions 

-0, 


(JL)±G° 


(3.27) 


to  remove  the  above  singularities  (see  Appendix  A) 
f  is  related  to  k  via: 


f  =  kh  if  z  in  v  =  khM,  M  ±  Di ag(  1  ,-a/c )  if  z  in  +  S4 

(3.28) 


(i.e.  fj  -  kjhj,  y2  =  k1h2 ,  f3 


k2hiM>  ^4  =  k2hlM^- 


Proof. 

(i)  We  first  note  that  the  product  of  the  jump  functions  g^  at  an  inter¬ 
section  point,  determines  the  nature  of  the  singularity  of  the  function 
y(z)  at  this  point.  For  the  sake  of  simplicity  assume  that  v  is  scalar 
a>id  that  it  behaves  like  zv  as  z  -  0,  in  Sj.  Then  y  <v,  z'^g^O),..., 

^4  x  zVg1(0)g2(0)g3(0) ,  this  implies  e~2lTTV  J  g(0).  Conversely,  if 

4  r  j  =  1 

n  g  (0)  =  e  ,  then  »  n,  zv  as  z  •*  0.  This  analysis  is  unique  within  the  trans- 
j  =  l  J 

formation  v  -*•  v  +  integer  and  zv  *  E  ^zvE,  and  its  generalization  to  the  case  that 
t  is  a  matrix  and/or  the  intersection  point  is  infinity  is  straightforward. 

In  what  follows  we  take  the  above  integer  to  be  zero  since  we  are  only  allowing 
f  to  have  an  integrable  singularity. 

Equation( 3 . 19)  implies 


j",  9J  '  1  j,  VM»  “  35  2  * 

and  hence  f(z)  %  Oiag(z  ,  z  )  (oQ  /  n/2).  Thus  y( z )  has 


a  monodromy  matrix  MQ  at  the  origin,  which  is  consistent  with  (3.21). 


-*  ao 


,  which  is  consistent  with  the  fact  that 


Similarly  n  g .  -v  M  as  ! z 
j  =  l  J 

'0*,  0^ 

f(z)  ^  Di ag ( z  ”,  z  °°)  as 


ao  # 


(ii)  Consider  the  following  transformations 


«1  =  kjhjAj. 


y2  =  ^ 1 ^ 2^2  ’  ^3  =  ^  2 ^ 2^3  * 


*4  '  k2hlA4 


(3.29} 


where  ,  j  =  are  constant,  non-s i ngul ar  matrices  and  the  functions 

k^,  k^,  h^,  h^  are  defined  in  +  S2 »  +  S^,  +  S^, 

respecti vely : 


Recall  that 
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Cr  kl  ■  V’lWi1*;1  fhjAjMg^-'h-1  on  C, 

or  ki  =  k2  | 

C3.  h^A^Mg^A^h,^  on 

(3.37) 

Letting  A^  =  A^  =  I,  =  1,  =  0,  equation  (3.37)  reduces  to  (3.23a) 

and  H  =  h.  We  note  that  the  k  problem  inherits  its  singularities  from 
the  f  problem:  Consider  the  product  of  the  jump  matrices  at  infinity  and 


2 


at  the  origin  (consider  (3.37)  with  =  A,,  =  I): 

9kJ9k1  ?  h292M  h2lhlMg4hl 1 ’  9k39k1  “  M-’ 

9k39kL  ^  hl(0)(hi1(0)h2(0))G2M'1(h21(0)hl(0))MG4M-hi1(0)’  Z  *  °- 


But , 


h^(0)h2(0)  = 


-a(p1(0)-o2(0>)  1 


since  Pj(0)  -  o2(0)  =  -1  (see  (3.23c)).  Also  M”*G^M  =  G^,  thus 

9'k^k1  ''  12I  *  *  •  9k|9ki'  hi(°)Eo‘MOlEOhl1(0)-  2  *  °- 

(3.38) 

Equation  (3.38)  implies  that  the  monodromy  matrix  of  K  at  the  origin  is 
h j ( 0 ) Eg E^h ~ 1 ( 0 ) .  This  is  consistent  with  the  facts  that  ^ h~ 1 

and  the  monodromy  matrix  of  is  Eq^M^Eq  (see  (3.14)). 

Equations  (3.38)  imply  that  k  has  the  same  singular  structure  as 
Y.  These  singularities  can  be  removed  by  using  the  auxiliary  functions 
(3.27). 


Proposition  3.3. 

Let  y(z)  be  the  solution  matrix  of  the  inverse  problem  formulated 
in  Theorem  3.1.  Then  y(t), 

y(t)  *  -(^  +  2t),  u  *  -21  im  Zf12(z)e'2q(z), 

|  z :  -*■<»> 

solves  PIV,  where  Yj2(z)  is  the  upper  left  entry  of  y(z). 


(3.39) 


Proof . 


Equation  (3.12b)  defines  u  in  terms  of  Y  (and  hence  in  terms  of  r), 
and  equation  (2.8b)  defines  y  in  terms  of  u. 


3.3.  Schlesinqer  Transformations. 

As  it  was  mentioned  before,  the  case  of  general  I'q,  ^  can  be  re¬ 
duced  to  the  case  of  0  <  Og  <  1 ,  0  <  0^  <  1.  In  this  section  we  pre¬ 
sent  the  transformations  which  shift  the  values  of  Og,  0^  by  half-integers. 
Similar  ideas  were  used  by  C.  Cosgrove  [13]. 


Proposition  3.4. 

Let  y  and  y'  be  solutions  of  PIV,  equation  (2.6)  with  a  =  20^-1, 
t?  =  -2Qg  and  a'  =  20^-1,  8'  =  -2(Og)2  respecti  vely .  Let  Y  and  Y‘  be  solutions 
of  the  corresponding  i somonodromic  problem  (2.7).  Consider  two  sets  of 
transformations : 


a : 


+ 


n 


G '  =  G  +  m 

CD  QD 


b: 


G0=  G0 


0-0 


2n  +  l 
2 

2m-*- 1 


m,n  e  2.  (3.40) 


Then : 

(i)  The  monodromy  data  for  Y  and  Y'  are  the  same. 

(ii)  The  solution  of  the  inverse  problem  for  Y'  can  be  obtained  from  Y: 
Y'  =  RY.  There  are  two  cases.  In  particular: 

R(z)  is  a  rational  function  of  z 
R(z)  is  z1/2 


(a) 

(b) 


times  a  rational  function  of  z. 


(3.41) 
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'  o  1  «  0  -  -I 

o  o  7* 


0 '  =  0  ♦  i- 

eo  ao  £ 


'e‘  •  on  *  l 

o  o  2 


O'  *  0 


1 

7 


y(v’V°» 

7  (v  -26q  ) 
(3.42a) 


(3.42b) 


(3.42d) 


The  transformati ons  (3.42)  generate  all  the  Schlesinger  transformations 
specified  by  (3.40).  For  example,  if  Rg,  Rg  are  defined  by 


R5 '  30  ~  0O’ 


OD  JO 


*  l’  R6:  =  '0  l’ 


Then 


R5  =  R1R3’  R6  ‘  R2R4 ' 


(3.43) 


The  above  transformations  naturally  induce  transformations  mapping 
solutions  of  P I V  to  solutions  of  P I V  with  different  values  of  the  para¬ 
meter.  For  example  R,  implies 


y1  =  M2T(t  -  ♦!],  u'  --  u ( “  -  *  - 


y 


•  =  .  2v  V)  .  (V  . 

y  y 


V 


-o 


:  3.441 
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where  Qq,  G^  are  defined  in  (3.42b)  in  terms  of  eQ, 

Proof. 

We  first  note  that  equation  (3.17)  is  invariant  under  the  transform¬ 
ations  Qrt  -  0'  Q  -  Q'  iff  Qi,  Q‘  are  given  by  equations  (3.40). 

(J  oo  od  U  00 

Then  Propositions  3.1,  3.2  imply  that  the  solvability  of  the  inverse 
problem  for  Y'  is  equivalent  to  that  of  Y.  Linder  the  transformations  a, 
b,  M  is  mapped  to  M1  where  M'  =  M  ,  M'  =  -M  respectively.  Thus, 

OO  OO  ODOD  OD  oo 

assuming  Y*  *  RY,  R  =  R.  if  z  in  S-,  then  the  jump  conditions  imply  a 

J  3 

RH  problem  for  R: 


a:  »jtl  .  Rj  on  Cj„  b:  Rj4,  *  Rj  on  CJt,  J  •  1,2,3. 

Rj  »  R4  on  Cj  Rj(z)  ■  -Ra(ze2'")  on  C,  (3.45) 

•  Equations  (3.45)  imply  (3.41);  to  determine  completely  the  form  of  R(z), 
i.e.  to  specify  the  rational  function  of  z,  one  uses  Y'  =  RY  to  obtain 
the  following  boundary  conditions  for  R: 

a:  R(z)  x  YQ(z)znI‘YQ1(z)  as  z  *  0,  n  ±  diag(l.-l) 

R ( 2 )  %  r (z)(i)mnY'1(z)  as  |z|  -  -  .  (3-46) 

b:  as  above  with  n  -  *  m  *  ^r~ ~  •  (3.47) 

If  y',  u',  v',  Qq  =  Gq  -1/2,  =  0^  +  1/2  are  the  transformed 

quantities  of  y,u,v,Qn,G  under  the  transformati on  given  by  R,(z),  i.e. 

u  °°  i 

Y  (z;y  ,u  ,v  tGg.G^)  =  R^ ( z ;y ,u , v  , Gq , O^) Y ( z ;y , . . . ) , 


1/2  are  the  transformed 


quantities  of  y\  u',  v'.  O'  V,  then  from  the  transformation  given  o.v 


i  .e. 


Y,,(z;y,,,u,\. . .  ,3")  =  R?(z;y‘,u' . 3x)Y‘  (z;y'  ,u' , 


a  tedious  but  straightforward  computation  shows  that, 


R~( z ;y ’ (y ,u , . . . ,0^) , . . . )R^( z ;y  ,u, . . .  ,Gx)  -  I 


Simi larly , 


R  (z;y' (y,u,. . . ,Gj ,. . . )R4(z;y,u 


=  I. 


A1  SO 


R^(z ;y' (y.u, ... ,0^) , )R4(z ;y ,u .... .o^)  R^(z), 


where 


Hence,  the  successive  application  of  the  Schlesinger  transformations 
defined  by  the  multiplier  matrices  Rk(z),  k  =  1,2, 3, 4,  maps  3q,  -x  to 
q'  =  Qn  +  n/2,  0'  =  3  +  m/2,  n,m  £  2.  To  obtain  equation  (3.44),  note 
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t ha t  the  above  transformations  map  equation  (2.7a)  as  follows 


y  =  AY  =  >  Y '  =  (RA  +  R  )R  Y' . 
z  z  z 


3.4.  Special  Solutions 

As  it  was  mentioned  earlier,  for  certain  choices  of  the  parameters 
a, 6,  PIV  admits  rational  solutions  or  one  parameter  family  of  solutions 
expressible  rationally  in  terms  of  the  Weber-Hermi te  functions.  Such 
solutions,  are  naturally  obtained  via  the  RH  formalism  presented  in  §3.2 
For  example. 

Example  3.1. 

Let  0  =  -Gn,  0  <  Q  <  i,  and  assume  that  a  =  c  =  0.  Then  b  =  - 

Eq  3  I,  and  the  solution  of  the  RH  defined  in  Theorem  3.1  is  given  by 

/d)°-  f 


'f(z)  = 


2  tt  i  j 


4  -  Z 


<±re. 


where  the  contour  C  is  defined  in  Figure  3.4  (C  is  below  the  branch  cut) 


Figure  3.4. 
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Associated  with  equation  (3.48),  the  solution  y  of  P I V  is  proportional 
to  He20  ^(it),  where  He^(t)  denotes  the  Weber-Hermite  function  [55] 

OD 

of  parameter  j. 


Proof 


Equation  (3.17)  implies  that  it  is  possible  to  choose  a  =  c  =  0 
provided  that  9^  =  *9q  +  n/2,  n  c  2.  Here  we  will  consider  the  case 
of  0  =  -en  and,  in  view  of  §3.3,  we  assume  0  <  '  <  1/2.  Equation 
(3.9)  implies  b  =  -d  and  =  diag(aQ ,a^ ) .  Using  the  similarity  argu¬ 
ment  (equation  (3.15))  we  take  a ^  =  1.  Thus  the  basic  RH  problem  re¬ 
duces  to 


cr 

3=1,  3  =  de2q 

,  =  1 

2 i " O  ,  .  2q  -2 1  - 

-2  i  -  ' 

(3.49) 

cr 

i=e  ®,  c  =  de  Me  r. 

r  =  e 

Letting  t  =  (^.t^)  the  above  reduces  to 


(3.50) 


The  solution  of  (3.50a)  satisfying  the  boundary  condition  (3.201,  is 
given  by 


Let  us  consider  the  homogeneous  problem  corresponding  to  (3.50), 


+ 


ri 


2  ' 


Its  solution  (satisfying  (3.20))  is 
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(3.52) 


Thus,  the  solution  of  (3.50)  is  given  by 


»<i>  ■  ♦uni  ♦jk  | 


(3.53) 


C3+Cl 


where  the  matrix  f  given  by  f  =  (0,  84'j)(<t  )  ,  i.e. 

fll  =  f21  =  f22  =  °’  =  d^-  &  + 


where  -  and  +  denote  the  limits  of  z  from  the  t  (S^-^)  and  -  (S-^+S^) 
regions  respectively.  Thus 

f12  =  de2q<^U)2j”  *'3*  f12  ■  de2<,(-^)2j’  Cr 

re  re 

,  Hence  substituting  the  above  in  (3.53)  and  using  (3.52)  we  obtain  (3. *18) 


4,  PAINLEVg  V 

The  fifth  Painleve  equation  (4.3)  can  be  obtained  [15]  as  the  compat¬ 
ibility  condition  of  the  following  linear  systems  of  the  equations 


Y  (z)  =  A ( z ) Y ( z ) , 


vt(z)  =  8; :  iy.'z; 


where  , 


A  ±  ♦  A„  -~r ,  B '  z  ' 

i  Z  C  Z  ”  * 


» 
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W  i  V  ♦  y(0Q  +  0^). 

The  compatibility  condition  of  (4.1)  gives: 

t  jj£-  ty  -  2v(y-l)2  -  ^(y-l)[(0o-8/0jy  -  (30o*0/0j]. 
t  SI  .  yv[,  ♦  ^0o-S,*0jJ  -  i(»*0o)[»  *  j(00*8,‘9.)]. 

t  SE  =  u{-2t-o0*y(v  *  7<VV0«^  *  *  jfOo*0,*8.)])’ 

(4-3) 


(4.2a) 

(4.2b) 

(4.2c) 


wi  th , 


The  general  fifth  Painleve  equation  with  non-zero  *  is  reduced  to  the 
case  6  =  -1/2  by  scaling;  the  case  6  =  0  may  be  transformed  to  the  third 
Painleve  equation  which  will  be  considered  elsewhere. 


ft 


I 


» 


ft 


i 


4.1.  The  Direct  Problem 
Propositi  on  4.1 


ft 


Let  Y(0)’  Y(l)  be  the  solutions  of  (4-la)  analytic  in  the  neighbor- 
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hood  of  2  =  0  and  z  =  1  respectively  and  normalized  by  the  requirements 
that  det  Y^  =  det  Y^j  =  1,  and  that  Y^,  Y^  also  solve  (4.1b). 

Let  Y^,  Y^  be  solutions  of  (4.1a)  analytic  in  the  neighborhood  of  in¬ 
finity  such  that  det  Y.  =  det  Y,  =  1  and  Y.  %  Y  as  |z|  - 
in  S j ,  where  Y^  is  the  formal  solution  matrix  of  (4.1a) 
in  the  neighborhood  of  infinity  and 

V  -  2  -  ar9  2  <  J.  S2:  i  -  ar9  2  <  T"  (4.5) 

The  contours  C^,  C^,  are  defined  by  arg  z  =  -  arg  z  =  arg  z  =  0 
and  0  <  Re  z  <  1  respectively 


Figure  4.1 


Then  the  analytic  functions  Y^,  Y(i)’  Yp  Y2  satisfy: 

(l)  Y(o)(z)  %  Y(0)(z)z  °  as  2  *  0.  %  •  Diag(^,-^),  GQ  f  n,  (4.6) 

where  Y^(z)  is  holomorphic  at  z  =  0  (Y^  has  a  logarith¬ 
mic  singularity  if  Oq  =  n). 
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G1  31. 


(ii)  Y(1)(2)  %  Y^(z)(z-1)  as  z  -  1,  Dj  *  Diag(  — — ),  f  n,  (4.7)  j 


where  Y^(z)  is  holomorphic  at  z  =  1  (Y^  has  a  logarith¬ 


mic  singularity  if  0^  =  n). 


(iii)  Yj(z)  ^  Yoo(z)e^z^)  “  as  |z|  ■+  -  in  S  ■ ,  Db  t  Diag(j^,-  ,  j  =  1,2, 

Q(z)  #  Diag(q.-q),  (4.8) 


X  t 

q  *  2~,  Y^z)  is  holomorphic  at  z  =  °°. 


2  i  -rr . 


(iv)  Y(0)(ze  )  =  Y(0)(z)Mg.  4 s  z  *  0 ,  Mq  f j 


17TOo  .  ’"0o 

e  2i  i  JQe 


-  i  1T&- 


(4.9) 


JQ  =  0  if  0Q  t  n,  JQ  =  1  if  Og  =  n. 


i  10, 


1  10, 


2ii, 


(v)  Y(i)(2e  )  =  Y(i)(z)Mi»  as  2  *  1*  Mi  *1 


2iiJ1e 


-  1  10, 


(4.10) 


J^=0  if  0j  f  n,  =  1  if  Oj  s  n. 


(vi)  Y  j  ( z )  =  Y^zjG^  Y^z)  =  Y2(ze2l'r)G2Mgo, 


(4.11) 


where 


1  0 


.  a  1 


g2  * 


1  b\  i iG^  -iiG. 

,  M  t  Diag(e  ”,  e 


(4. 12) 


■v11^  Yi  =  Y(o)Eo’  det  Eo  "  1;  Yi  =  Y(i)Er  det  Ei  =  1- 


(4.13) 


)  v  i  i  i )  Let  Y^,  Y*  dentoe  the  limits  of  Y^  as  it  approaches  from  above 


and  below  respecti vely .  Then 


’I  ■  yTe11miei 


(4.14) 


I 
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I 

I 

I 


Furthermore,  the  parameters 

MO  i  {a,b,a0,  8q  ,  yq,  <Sq,  o1.  Bp  Yp  ip', 

i.e.  the  entries  of  Gp  Gj,  Eq,  Ep  satisfy  the  following  con¬ 
sistency  conditon. 


(ix) 


Proof. 


w  •  E3VEoE;V£r 


(4.15) 


1.  Analysis  near  2=0: 

Let  Y(0)  =  ^Y(0)’  Y(0)^’  then  ^or  <  1  we  find: 


y(l),,,  .  y2  »  /°0, 
r(0)<2>  -  z  < 


/o\  *8  / 2  0 

r<Ej(Z)  -  2  0  e  0 


Z  +  . 


i ,  Q0  t  n 

(4.16a)) 

(4.16b) 


where . 


&  •■>('•  y 

■y*^- A- ?>•«*>«• 

%(')  ♦  f  t-jk « - 1<  w ♦ 


Expressions  for  Kq\  may  be  given,  but  they  are  not  needed  in  the  discussion.  i_ 

The  constants  with  respect  to  z  multiplying  the  {  •  above  are  fixed  by 

the  requirement  that  equations  (4.16)  satisfy  (4.1b).  We  note  that  when 


i 


Oq  =  n  there  will  be,  in  general,  a  logarithmic  term  and  the  two 
linearly  independent  solutions  are 


Y(0),2)-  V(0)(2>  ■  Jo(tn  2),(0)(2>  *  2  °°  V(0)(2)’  ■0  ■  n  <4-17> 

-  (2) 

where  is  a  complex  constant  and  Yi  (z)  is  a  polynomial  in  z.  Equa¬ 


tions  (4.16)  imply  Y^(ze21lT)  =  YQ(z)e 
(4.16),  (4.17)  imply  (4.9). 


2  i  -jt  Dr 


Simi larly  equations 


2.  Analysis  near  z  =  1: 


Let  Y^j  =  (y|||,  y|2|),  then  Yor  i z- 1 1  <  1  we  find 


,(1) 


Y  (Z)  =  (z-1) 


(1) 


(1) 


(z-1)*  ... 


(4.18b) 


where , 


6 


rJ<“ *  r><^>  *  ; 


(2)  Hw-e  /2  w+6  /2  0  l+2w 

Ll  *  “1-6/  [y (T+w^0  J’(v  +Go )  ♦  vy(?  *  V  ♦  r)( TTw^7l)3. 


0,(t)  ♦ 


rt  ■■  1,..  .  61n  1 


(fr{v  -  -(w  ♦  j-)]  -  7 


)df 


Y(l)(z)-  Y(l)(z)  =  Ji[Mz_1)]Y(l)(z)  *  (2-1)  "l/2Yjij(z),  =  n, 


1 


(4.19) 


-  f  2 ) 

where  is  a  complex  constant  and  Y|^(z)  is  a  polynomial  of  z-1. 
Expressions  for  may  be  given  but  they  are  not  necessary  in  thls 


discussion. 
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3.  Analysis  near  2  = 

The  two  linearly  independent  formal  solutins  of  (4.19)  have  the 
expansions 


(2) 


Y  ( 2 ) 


(z) 


where , 


(4.20a) 


(4.20b) 


K  ♦ 


it* 


-L(w  ♦ 

y' 


^)][v 


♦  6. 


y(w  ■  2^-)]  -  w , 


Q 


Using  similar  arguments  to  those  used  in  Proposition  3.1  we  obtain 
(4.11). 

Let  Jj  denote  the  monodromy  matrix  of  Y^  at  z  =  1.  Then  Y*  =  Y~J^. 
However,  =  Y^E^  and  is  the  monodromy  matrix  of  Y,^.  Hence 
Jl  =  which  implies  (4.14). 

(4)  Consistency. 

Let  and  Jg  denote  the  monodromy  matrices  of  Y^  at  z  =  1  and  z  =  0 
respectively.  Equations  (4.11)  near  z  =  0  imply: 


Y'(z)  =  Y2(ze2in)G2M^ 

Thus  G^M,  = 

(4.15)  follows. 


Y^(ze  )G^G2Mao  =  Y  ^  (  z )  JgG^G2M;r  Y  ^  (  z  )  J  j  JgGjG^^. 
But  Jj  =  E^M^Ep  JQ  1  E^MqEq,  thus  equation 


I 


I 
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Proposition  4.2. 

(i)  The  monodromy  data,  MD,  defined  in  Proposition  4.1  are  time  in¬ 
variant. 

(ii)  All  of  the  MO  can  be  expressed  in  terms  of  two  of  them.  This 
follows  from: 

1.  det  Eq  =  det  E^  =  1. 

2.  Equation  (4.15). 

3.  If  Y  solves  (4.1)  with  y  satisfying  PV,  then  Y  t  R'^YR, 

1/2  -i /? 

R  t  Diag(r  ,r  ),  r  constant,  also  solves  (4.1)  with  y  satisfy- 
in  PV.  The  Stokes  matrices  and  the  connection  matrices  Eq,  E^ 
are  transformed  to  Gj  =  R~ 1 G ^ R ,  E-  =  R'^E.R,  j  =  1,2,  i  =  0,1, 
i  .e. 

a  =  ra,  b  =  b/r,  a.  =  a.,  s.  =  i./r,  y.  =  ry  .  ,  6.  =  i.,  i  =  0,1.  (4.21) 


4. 


Changing  the  arbitrary  integration  constants  :n,  :, 


amounts  to  multiplying  y|||(z),  y|^|(z)  by  p.  and  p*^,  i  =  0,1, 

A  1 

respectively.  This  maps  E-  to  E.  =  P.E^ ,  P.  f  Diag(pi,  pT‘),  i.e. 


«i  *  Piai ’  •«  •  pi8i- *  r-  ‘i  *  r-  '  *  0>1- 


Proof. 

Similar  to  that  of  Proposition  3.2. 

4.2.  The  Inverse  Problems. 

In  what  follows  we  formulate  a  RH  problem  for  the  case  that 
0  <  0q<  2,  0  <  ©1  <  2,  0  <  <  2.  This  assumption  leads  to  a  regu¬ 

lar  RH  problem. 
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Theorem  4,1. 

Consider  the  following  matrix,  regular,  homogeneous  RH  problem 
along  C^,  C2,  (Figure  4.1):  Determine  the  sectional ly  holomorphic 
function  f(z),  f(z)  =  fj(z)  if  z  in  Sj ,  j  =  1,2  from  the  following  con¬ 
ditions 

1.  T.  satisfy  the  jump  conditions 

f2(c)  =  f1(c)g1(c).  YjU)  =  ’  fi^  =  (4-23) 

along  the  rays  C^,  C, ,  respecti vely ,  where 


h  s  «2  s  93  *  eQ*l\ he'Q-  (4'24> 

2.  f(z)  ^  (|)°-( I  ♦  0(£))  as  i z |  -  -  .  (4.25) 

3.  f(z)  has  at  most  mtegrable  singularities  at  z  =  0,  z  =  1,  and  z  =  * 
the  monodromy  matrices  of  ^  are  given  by 


-1  -1  ”^1 
E0  M0E0’  e  E1  MlEle  ’  M~’ 


;  4 . 26 i 


respeci tvely ,  where  Q.  =  Q(l).  In  the  above  6 ■ ,  Q,  M  ,  D  ,  Mn,  M.  are 

i  j  ®  30  u  i 

defined  in  Proposition  4.1. 


4.  The  monodromy  data  satisfy  the  properties  of  proposition  4 . 2 ( i i ) . 
Then : 

(i)  The  above  RH  problem  is  discontinuous  at  z  =  0,  1,  «.  Actually 


h  '  ^h*’  2'0;  9192 '  "*• 2 


-1  Q1P- 1M- lr  '°1 
g3  --  e  E,  M,  E,e 


1  '  i  "r 


4  _  2’ ) 


and 


x>  ' 


1. 
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(ii)  To  obtain  the  solution  of  the  above  RH  problem  consider  the  follow 
ing  RH  problem  along  the  contour  Cq  defined  by  Im  z  =  0: 

i  Im  2 

i 

i 


Figure  4.2 


Determine  the  sectionally  holomorphic  function  k(z),  k(z)  =  k ^ ( z )  if 
Im  z  >  0,  k(z)  =  k^(z)  if  Im  z  <  0,  from  the  following  conditions: 


1.  k^,  k^  satisfy  the  jump  condition 


h^Mh^ 


1 


C0:  kl  ‘  k2  ' 


h(  z : 


1 

vao(z) 


h.MM  g0h" 
1  1 

h.MM  hT1 
I  00  1 

0N 

1 


1 


(2)  = 


Re  z  <  0, 

0  <  Re  z  <  1  M  t ! 
Re  z  >  1 

I  f  dee 


-a/b 


2-ri  ! 


;-2  ’ 


Cl+C2 


and  hj,  h^  denote  h  for  Re  z  >  0,  Re  z  <  0  respecti vely. 

i  D  i 

2-  k(z)  ^  (j)  ”( I  +  0(~) )  as  !z  j  -  ®  . 


1  on  C,  ,  -1 

(4.28) 

(4.29) 


on 


3. 


k(z)  has  at  most  integrable  singularities  at  z  =  0,  z  =  1  with 


monodromy  matrices  given  by 

hl(°)E"1M0Eohi1(0) »  h1(l)eQ{1)E'1M1E1e'Q(1)h^1(l)  (4.30) 


respectively. 

The  above  RH  problem  is  discontinuous  at  z  =  0,  1,  °°.  Actually  if 

g.  ,  g.  ,  g.  denote  the  jump  matrices  for  Re  z  <  0,  0  <  Re  z  <  1, 

1  2  3 

Re  z  >  1  respectively  then 

gj^Xgk  n*  as  |z|  -  gk  1gk  ^  h^OE^M'^Qh'^O)  as  z  -  0; 

(4.31) 


g^gk  x  h1(l)eQ(1)Ej1Mj1E1e”Q(1)hJ1(l)  as  z  -  1. 
2  3 


However,  the  above  RH  problem  can  be  mapped  to  a  continuous  one  using 
the  appropriate  auxiliary  functions  (see  Appendix  B). 

Y  is  related  to  k  via 


Y*  =  klhr  ^ I  =  kihlMM®’  f2  =  klh2'  (4.32) 

Proof 

(i)  The  products  of  the  jump  matrices  at  a  given  point  determine  the 
nature  of  the  singularity  at  this  point.  Equations  (4.27),  which  follow 

V2  -'0/2 

from  equations  (4.24),  imply:  y(z)  ^  D i a g ( z  ,  z  ),  z  -  0,  Qn  /  n 

°1/2  "V2 

and  Y  has  also  a  log  z  term  if  ©q  =  n;  y(z)  ^  D i a g ( ( z - 1 )  ,  (z-l) 

and  Y  has  also  log  (z-l)  term  in  =  n;  y(z)  ^  Diag((-)'®,  (-)  as 

|  z  |  -  ®.  Hence  y(z)  is  singular  at  z  =  0,  1,  .  .  Assuming  0  ~q  1, 

0  <  Q.  <  1,  0  ^  G  <1,  y(z)  is  integrable  at  the  above  points. 


Figure  4.3 


(ii)  Consider  the  following  transformations 


'l  =  klhlAr  y  1  =  k2hlA2  ’  f2  =  k lh2A3  ’  f2  =  k 2 h 2 A4 ’ 


(4.33) 


'  where  is  if  1  y  0  and  if  Im  z  <0  (clearly  ^  in 


however  we  use  this  artificial  separation  in  order  to  make  the  h-RH 


problem  continuous) 


Recall  C2 ■  ^2  f|9i*  Y^  -  ^3’  ^1  ^1^3 


(4.34) 


Equations  (4.33),  (4.34)  imply: 


C2:  h2  =  hlAl9lA3  ’  Cl:  hl  h2A4g2A2 


We  choose  the  A's  in  such  a  way  that  the  h-RH  problem  is  continuous  both 


at  zero  and  infinity. 


Continuity  at  zero:  I  as  z  -  0, 


Continuity  at  infinity:  A^g jA^ ^A^A,, *  -  I  as  z 
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or 


or 


WzWJ?  -  '■  VIVa' 


M  *  A^A,  =  A^A.M'1,  G,  MG-  =  M,  i.e.  M 
4  4  14®  14 


0  M 


12 


(4.35) 


b  M12  M22' 


(assuming  a.b.M^  /  0).  Hence  the  h-RH  problem  becomes 


h2  -  hj 


A1g1A3  on  C2 
LA.MM^g^M^A"1  on  Cj 


,  1 

or  =  hjAj ^  ^Axt 


A-i 

3  ’ 


(4.36) 


where  X  =  -1  on  C0,  X  =  1  on  C^.  Let  H,  =  h0A,,  H,  =  h1A1  and  (4.36) 


2  "23’  1  11 


reduces  to 


C!  +  C2 '  H2  =  Hl\ 


1  0 


xxt  , 

.ae  1 


1  on  Cp  -1  on  C^. 

(4.37) 


Since  the  H-RH  problem  is  continuous  at  ®  we  look  for  a  solution  such  that 

H  %  I  as  z  -  ®. 

(H^,H^2))  =  (H^1 )  ,h[2)  )  +  aem(H|2),0), 


or 


Thus 


H<2>.(°),  h<»  -  H<»  .  ae'^(°). 


1  0 


,ap  1 


,  p  as  in  (4.28) . 


(4.38) 


^1  =  ^1  on  2  ax’s»  for  1  '■  however  equations  (4.33a) 
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(4.33b)  imply  that  k^  f  k ^  for  Re  z  >  1  (since  A^  f  A^ ) .  Hence,  al¬ 
though  has  a  discontinuity  only  along  C ^  (0  <  Re  z  -  1),  k  has  a 
discontinuity  along  the  entire  positive  Re  z  axis.  Similarly,  equations 
(4.33c),  (4.33d)  imply  that,  although  ^  is  continuous  m  $„,  k  ’s  <j’$- 
continuous  in  S^;  we  choose  the  negative  Re  z  axis  to  be  the  c.r.e 
the  discontinuity. 

Let  us  now  formulate  the  K-RH  problem: 


IM  z  =  G:  ^2  =  * i'  z  '  0;  ^  =  ^g7,  0  ■  Re  z 


Using  (4.32)  we  find 


Im  z  =  0:  k.  =  k. 

1  £ 


Using  A^  =  A^M,  A^ 


kl  =  k2  < 


^2^4^3  ^2 

Re  z  ••  0 

h 1 A2g3Al 1 h  i 1 

0  •  Re  z  •  1  . 

4  .  39 

Re  z  •  1 

we  find 

Re  z  <  0 

0  <  Re  z  <  1 

.  h, A,MM  AT 1h7  L 
\11  00  1  1 


Re  z  >  1 


which,  with  the  choice  of  A^  =  A^  =  I  ,  =  1,  =  0,  reduces  to 

(4.28).  The  k  problem  inherits  its  singularities  from  the  *  problem. 
Consider  the  product  of  the  jump  matrices  at  the  singular  points: 


g^  .  h^Mh^  ,  g^  •  3  "  ' 


Skjgfc  as  '2 


9kjsK2--  hl'0»Eo’,Ho:EoS1(0>  «  2  •  0, 


Then 
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-  h1U)eqC1)E-1M-!E1e-Q(1)h-;(l) 


2  3 


as  z 


(4.401 


Equations  (4.40a),  (4.40c)  are  obvious.  To  obtain  equation  (4.40b) 
note: 


hT1(0)h?(0)  =  ( 

1  2  \a(p2(0)-Ol(0)) 


=  G1 ,  since  o2(0)-c1(0)  =  1 


( see  (4.38)). 


Hence , 


q‘1gk  =  h2(0)M'1h'1(0)h1(0)MMoog3(0)hJ1(0)  =  h  3  ( 0 )  G  1M_  1 1  MM^g  3  ( 0 )  h  * 1  ( f 
1  2 

=  h1(0)G1G2MaoE‘1M1E1h‘1(0)  =  ^  ( 0 )  E*  1  ( 0) ,  (4.41) 

where  we  have  used  M  ^G^M  =  G2  (see  (4.35)),  and  (4.15). 

Equation  (4.40b)  implies  that  the  monodromy  matrix  of  k  at  the 
origin  is  h  ( 0 ) E^M^h" 1  ( 0 ) .  This  is  consistent  with  the  facts  that 
k  =  h j  and  the  monodromy  matrix  of  fj  is  Eq  MqEq-  Similarly  for  the 
monodromy  matrix  of  k  at  z  =  1.  Thus  k  has  the  same  singularities  as  f. 
These  singularities  can  be  removed  by  using  appropriate  auxiliary 
functions . 


Proposition  4.3. 

Let  f(z)  be  the  solution  matrix  of  the  inverse  problem  formulated 
in  Theorem  4.1.  Then  y(t),  which  may  be  obtained  from  (4.1a),  (4.2a), 


solves  PV. 


4,3.  Schlessinqer  Transformations. 

Proposition  4.4. 

Let  y  and  y'  be  solutions  of  PV  equation  (4.4)  with  i,  2,  and 
a',  S',  y'  respectively,  where  a,  3,  y  are  related  to  Oq,  via 

(4.4).  Let  Y,  Y*  be  solutions  of  the  corresponding  i somonodrorm c  prob¬ 
lem  (4.1).  Consider  the  sets  of  transformati ons : 


Q0  =  G0  1  n 

r>,  •  - 

"0 

30 

-  i 

'0 

-o  -  n 

a : 

9i  =  °1 

b:  0'  = 

H  i  n 

c : 

-  r 

1 

=  Oj  +  m, 

O'  =  G  +  m 

GO  OO  - 

-  i  _ 

~  OO 

+  m 

ac  — 

OO 

OO 

where,  n,m  are  either  even  or  odd  integers.  (4.42) 

Then : 

(i)  The  monodromy  data  for  Y  and  Y‘  are  the  same. 

(ii)  The  solution  of  the  inverse  problem  for  Y'  can  be  obtained  from  Y: 

1/2 

Y'  =  RY;  a:  R(z)  is  z  times  a  rational  function  of  z, 

1/2 

b:  R ( z )  is  (z-1)  times  a  rational  function  of  z,  (4.43; 
- 1/2  1/2 

c:  R(z)  is  z  (z-1)  times  a  rational  function  of  z. 

In  partiuclar: 


( 


1 


o: 


V 


0' 


Q  +  1 


[j<*  * 


-1,  2 

z 


(4.44a 
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The  transformations  (4.44)  generate  all  the  Schlesinger  transformations 
specified  by  (4.42).  For  example, 

Rj  +  ^Rj  =  I>  j  -  1 , 3 , 5 , 7 ;  ^1^7  ~  »  ^2^8  ”  ^10’ 

where  • 

Rg:  s i  ■  0„  ♦  1,  e;  =  0;  *  l,  o:  ■  0.; 

R10:  0O  s  0O  Qi  =  V1,  =  0®‘  (4-45) 

The  above  transformations  induce  transformations  mapping  solu¬ 
tions  of  PV  to  solutions  of  PV  with  parameters  related  as  in  equations 
(4.42). 

Proof 

Equation  (4.15)  is  invariant  under  the  transformations  3q  -  0^ , 

Ql  -  0^,  -*•  0^  iff  the  0's  transform  as  in  (4.42).  It  will  turn 

out  that  it  is  sufficient  to  consider  m  =  n  =  1.  Since  M, ,  M  -  *  M.  , 

i  °°  i 

assuming  Y'  =  RY,  equations  (4.23)  imply 


a : 

R2  R 1 » 

B:  R2  »  R*  , 

c : 

R2  =  R1 

on 

C2 

»; .  r- 

r1  =  -r; 

ri  =  -r; 

on 

C3 

Rj  =  -R2 

R1  =  ”R2 

R1  =  R2 

on 

C1 

I 


Figure  4.4 


Equations  (4.46)  imply  (4.43);  to  determine  completely  the  form  of 
R(z),  use  Y1  =  RY  and  the  boundary  conditions  for  Y  to  obtain  appropri¬ 
ate  boundary  conditions  for  R  (the  details  are  given  in  [45]). 


3.4.  Special  Solutions 


Example  4.1. 


Let  ©„=  -(0Q  +  ©j),  0  <  0O  <  1,  0 


<  s,  <  1  and  assume  that 


a  =  0.  Then  6^  =  Yj  =  Yq  =  0,  and  the  solution  of  the  RH  defined  in 


Theorem  4.1  is  given  by 


Hz)  - 


O  J2  0,/2 

z  0  (z-1)  1 


bz'0°/2(2-i)-&l/2_L  [  dcectc  °( ;-l )  1 
DZ  lZ  i]  2tt i  Jc  ;  -  z 


-3/2  -0/2 

Z  0  (Z-1)  1 


(4.47) 


where  the  contour  C  is  along  in  the  region.  Hence  associated  with  eauatii 
(4.47),  is  a  solution  y  of  PV  which  is  proportional  to  W  .  ,  ,,  where  W 


denotes  the  Whittaker  function  [55] . 


:'0  ’  ^  1  ( 1 ) 


To  derive  the  above,  note  that  in  this  case  the  basic  RH  problem 


reduces  to 


Figure  4.5 


+ 

Y  =  f 


Im  z  =  0 ,  Re  z  '  1  :  i=  1 ,  y  =  1 ,  r  =  0 

-1"-!  ”  i 

’  Im  z  =  0 ,  0  '  Re  2  <  1 :  a=  e  ,  v=  e 

-  iiO  _  -i-G 

C,  :  a=  e  ®,  y-  e  *. 


+;t 


Letting  f  =  (fj,  f.,)  the  above  reduces  to 


(4.48) 


f 2  =  ,  *0  =  y^  +  Bfj 


(4.49) 


Equaiton  (4.49a)  implies 


qq/2  S^/2  i 

2  U  (Z-1)  1  (J) 


(4.50) 


The  "homogeneous"  version  of  (4,49)  yields 


,  V2  V2 

2  0  (2-1)  1 


-0n/2  -0,/2 
2  0  (2-1)  1 


Thus  the  solution  of  (4.49)  is  given  by 


(4.51) 


y(z)  -  *(Z)(I  *  ^  j 

L1 


(4.52) 


where  f  =  (0,  Bfj)(<i>  )  ,  i.e. 

-i^G+U  V2  :■  n  /  2  :J  2  ,/2 

fll  =  f  12  =  f22  =  °’  f21  =  be  [4  U  ( 4 *  1 )  1  (4-1)  1  ]  + 


where  -  and  +  denote  the  limits  of  z  in  and  respectively.  Since 

r t  ®0^2  0./2  2 

f 2 1  =  u  U  (4-1)  1  )_  we  obtain  (4.47). 


APPENDIX  A 


In  Theorem  3 . 1  we  mapped  the  basic  p*  pod  i°m  which  under' 


les  the  initial  ,  alue 


problem  of  PIV  to  a  simpler  RH  problem,  i.e.  equations  ( 3 . 23  )- ( 3 . 26 ) 


'his  -  e 


problem  is  discontinuous  both  at  zero  ano  ’ne 
of  the  jump  matrices  at  z  =  0  and  z  -  is  given  by 


actual ly  the  product 


:  A .  l ) 


respectively.  We  now  map  this  discontinuous  problem  to  a  continuous  one. 
The  basic  idea  is  to  use  appropriate  discontinuous  auxiliary  functions 
such  that  the  product  of  the  jump  matrices  of  the  transformed  problem  is 
I.  This  procedure  is  the  well  known  [46]  so  the  derivation  is  omitted 
(details  are  given  in  [45]). 

Proposition  A.l. 

% 

Consider  the  k-RH  problem  formulated  in  Theorem  3.1  and  defined  by 
equations  ( 3 . 23 ) - ( 3 . 26 ) .  Assume  0  <  %  -  1,  0  <  ;  -  1,  -  /  \/2 

Define  the  sectionally  holomorphic  function  4>(z),  t  ( z )  =  }*(z)  if  z  is  in 
Sj  +  S2>  «(z)  =  $+(z)  if  z  is  in  S3  +  S4  as  follows: 


kj(x),  8j(x) 

Xq(x),  X”(x) 


k,(x).  §*(x) 

xo(x).  XJx) 


Figure  A. 1 


k^z)  =  ♦  (z)Xoo(z)RooXq(z)R0,  z  in  Sj  +  S., 


(A. 2) 


k2(z)  =  *  +  {z)X*(z)RjJ(z)rJ,  z  ^  S3  +  S4, 


where, 


z  ,°0  ,  z  v_0Ot  v+ 


■  DiagC^)  (jTy)  UL  xj{z)  ■  Oiag[(jfy)  °,  (j^f'0],  (A. 3) 


1  /e°°l  „  + 


X‘(z)  =  Oiag[(~)  ,  (yfy)  ],  Xjz)  =  DiagC ( jry)  ,  (jjj)  "].  (A. 4) 


The  auxiliary  functions  Xq,  X^  are  defined  with  respect  to  a  finite  branch 
cut  from  z  =  -1  to  z  =  1,  and  an  infinite  branch  cut  from  z  =  1  to  z  =  -1 
respectively;  these  branches  are  specified  by 


1  im  X‘(z)  =  I ,  1 im  X’(z)  =  I . 

Izl  —  0  z-0 


(A. 5) 


The  constant  matrices  RQ,  R^  are  defined  by: 

Rq  ?  E0h^(0),  Rj  t  E0G1G2G3M~1h’1  (0) ,  R~  f  (Rq)'1,  R*  *  (RqM)"1  ,  M  =  Diagd 


(A. 6) 


Then  the  <*> - RH  problem, 


VC3:  «+U)  =  4>'(c)x’1(£;)Rlx-(c)R'gk(r )[x^(c)R>o( c)Rq]_1  ,  (A. 7) 


4>’(z)  -  I  as  |zj  -  <»  in  S1  +  S^, 


where  is  the  jump  matrix  of  the  k-RH  problem  given  by  (3.23),  is  con- 
ti  nuous . 


Remark  A. 1 . 


The  case  Gn  =  0  can  be  handled  in  a  similar  way:  In  this  case  k 


1  +  i 


Xq,  X1  are  defined  with  respect  to  finite  branch  cuts  between  z  -  Zq  * 
and  z  =  Zq  =  1  -  i  passing  through  z  =  0  and  z  =  1  respectively;  is 
defined  with  respect  to  an  infinite  branch  cut  between  the  points  z  =  Zq 
and  z  =  Zq.  These  branches  are  normalized  by 


lim  xqU)  =  I*  lim  Xj(z)  =  I,  lim  Xa(z)  =  q,  lim  Xoo(z)  =  X^ 
|z|+“  |z|-*-°°  z-*-0  '  ’  z-*-l  ’ 

(B.5) 

with  [XJ-O^XJ-)  =  M^.  If  the  constant  matrices  are  appropriately 
chosen  then  <t> ( z )  satisfies  a  continuous  RH  problem  (the  details  can  be 
found  in  [45]). 


Remark  B.l 

The  cases  ©q  =  0,  =  0  can  be  handled  in  a  similar  way.  The 

logarithmic  singularities  can  be  removed  by  using 

'l  £n  F' 

1 

where  Zq,  Zq,  z^,  z^  corresponds  to  F 

(— — ).  (7^7-h  (f -7-)  respectively. 

z-i0  z'z0  z-iQ  z’z0 


X  ( z )  = 
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